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Context. Measurements of the 2.7 K cosmic microwave background (CMB) radiation now provide the most stringent constraints 
on cosmological models. The power spectra of the temperature anisotropics and the ii-mode polarization of the CMB are explained 
well by the inflationary paradigm. The next generation of CMB experiments aim at providing the most direct evidence for inflation 
through the detection of S-modes in the CMB polarization, presumed to have been caused by gravitational waves generated during 
the inflationary epoch around lO^^^'s. The S-mode polarization signals are very small (<10"^K) compared with the temperature 
anisotropics (~ 10"*K). Systematic effects in CMB telescopes can cause leakage from temperature anisotropy into polarization. 
Bolometric interferometry (BI) is a novel approach to measuring this small signal with lower leakage. 

Aims. If BI can be made to work over wide bandwidth (~ 20 - 30%) it can provide similar sensitivity to imagers. Subdividing the 
frequency passband of a Fizeau interferometer would mitigate the problem of 'fringe smearing.' Furthermore, the approach should 
allow simultaneous measurements in image space and visibility space. 

Methods. For subdividing the frequency passsband ('sub-band splitting' henceforth), we write an expression for the output from every 
baseline at every detector in the focal plane as a sum of visibilities in different frequency sub-bands. For operating the interferometer 
simultaneously as an imager, we write the output as two integrals over the sky and the focal plane, with all the phase differences 
accounted for. 

Results. The sub-band splitting method described here is general and can be applied to broad-band Fizeau interferometers across the 
electromagnetic spectrum. Applications to CMB measurements and to long-baseline optical interferometry are promising. 

Key words, cosmology: cosmic background radiation — techniques: interferometric — instrumentation: interferometers 
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I [ Interferometry has a long history for astronomical measurements 
. at radio, millimeter, submillimeter, IR and optical wavelengths. 
\^ In 1868 Hyppolyte Fizeau (Fizeau 1868) described how the di- 
ameters of stars could be measured by optical interferometry. 
Q^ He proposed 'masking' the aperture of a telescope to create in- 
lO terference fringes in the focal plane, similar to a Young's double 
\^ slit interference experiment. Today aperture masks are used to 
overcome atmospheric 'seeing' effects to reach the diffraction 
limit of single aperture telescopes (Tuthill et al.i2000 ). For even 
I , greater angular resolution, beam combination from widely sep- 
ILJ ' arated apertures is used in long baseline optical interferometers. 
. !^ ' Michelson used this technique to measure, for the first time, the 
diameters of stars. He used flat mirrors to reflect the beams from 
$_i separated apertures into a single telescope, which acted as the 
5^ beam combiner. 

This type of beam combination, in which beams are com- 
bined in the image plane of a telescope, is called 'Fizeau inter- 
ferometry,' or 'adding interferometry.' An alternate approach in- 
volves combining the beams in the pupil plane of a telescope and 
is called 'Michelson interferometry' after the technique used in 
the Michelson-Morley experiment (Monnier 2003). Both types 
of combiners are used in long-baseline optical interferometers 
(iTraubl 120001) . Conventional radio interferometers can also be 
thought of as Michelson interferometers. They typically mix the 
RF signals from each antenna in an array to lower frequencies 
(heterodyning) and interfere (multiply) them electronically one 
pair at a time with either an analog or digital correlator to mea- 
sure visibilities. This approach is sometimes called 'multiplying' 



interferometry. These techniques are widely used for radio wave- 
lengths to sub-mm wavelengths. As the number of antennas, A^, 
increases the number of correlations to be performed grows as 
A^(A^ - l)/2. Although radio correlators are improving rapidly, 
they are currently limited to combining signals from about 10 
antennas and bandwidths of a few GHz (iLawrence et al ] l2008h . 

In this paper we describe a Fizeau 'adding' interferometer 
that overcomes the bandwidth and large-A^ limitations of hetero- 
dyne interferometers. The instrument is optimized for precision 
measurements of the temperature and polarization anisotropy 
of the CMB a nd is sometimes called a 'bolometric interferom- 
eter' (Timbie et all 1200 6; Tucker etal. 2008 
2008; 



■Hamilton et all i2008. 



Charlassier et dj 
■Hamil ton & Charlassie^ ,201(il) . 
Interferometers are less sensitive to some kinds of systematic ef- 
fects found in imaging instruments (lBunnll2007i) . The technique 
can be used at any wavelength. 

In particular, we focus on an approach to broadening the 
bandwidth of a Fizeau interferometer. Spectral resolution is not 
required for many applications where the source has a contin- 
uum spectrum (such as the CMB) and signal averaging over 
broad bands is required to uncover faint signals. As one increases 
the bandwidth the different wavelengths within the passband will 
smear out the interference fringes and reduce the sensitivity of 
the instrument. Furthermore, this 'fringe washing' reduces the 
resolution of the interferometer in the u-v plane. In this paper 
we present a simple and powerful technique, which we call 'sub- 
band splitting' to overcome the limitations caused by large band- 
widths in Fizeau interferometers. We discuss how this approach 
can improve parameter estimation in the specific case of obser- 
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vations of the CMB, where sensitivity and u-v space resolution 
are critical for constraining CMB power spectra. 

We begin by motivating the need for adding interferometry 
in ^ followed by the reason sub-band splitting is required. ^ 
develops the formalism to describe how visibilities are measured 
by a Fizeau combiner (AppendixlAlpoints out how this arrange- 
ment can be used as an imager), ^describes in detail the tech- 
nique for solving for visibilities in sub-bands. 



2. Cosmological Context: The CMB and Adding 
Interferometry 

We now have a 'standard model' of cosmology in which the in- 
flationary paradigm describes many aspects of the universe ac- 
curately through well-quantified parameters. Precision measure- 
ments of the CMB are the most powerful probes for determin- 
ing these parameters. In particular, measurements of the angular 
power spectrum of the CMB temperature and E-mode polariza- 
tion anisotropy have yielded a wealth of information about the 
early universe (see Komatsu et al. 2009, for example). However, 
inflation is driv en by physics that we do not currently understand 
dBaumann et al.„2009.) . The most direct way to probe inflation is 
through the so-called fi-modes in CMB polarization. The am- 
plitude of the B-mode signal is directly related to the energy 
scale of the particle interactions that occurred during inflation. 
However, t he energy scale of i nflation is not know n. Recent mea- 
surements (iGupta et al.ll2010HChiang et alJl2010 ). are approach- 
ing the level at which B-modes are expected to appear in models 
in which inflation occurred at the GUT scale. However, inflation 
may have involved lower-energy interactions; there is no lower 
bound to the amplitude of the g-mo des. 

Current dTakahashi etap 120101: iHinderks et all l2009 l) and 
planned CMB instruments ( Bock et aP 2009h all use some type 
of imaging technique. With focal-plane arrays of hundreds of 
background-limited detectors they are capable of detecting the 
fi-mode signals predicted in the most optimistic models, at the 
level of ~ 10"** K. Systematic effects have been extensively 
studied for imaging polarimeters in the context of CMB mea- 
surements and appear to be con trollable at this level as well 
dHu et al.ll20"03HBock et"al]|2006h . However, at some level all in- 
struments can 'mix' the relatively large temperature anisotropy 
and E-mode polarization signals into Z?-modes. 

Different systematic effects are found in interferometers. It 
is for this reason that heterodyne interferometers have been used 
for many years to study the CMB temperature and polarization 
power spectra and the Sunyaev-Zel'dovich effect. These instru- 
ments multiply together the RF signals from all possible pairs 
of antennas (baselines) in an array to measure a set of visibili- 
ties, which are related to the sky image through a Fourier trans- 
form dRohlfs & Wilsonl2004l ). In fact, the first detection of CMB 
E-mode p olarization was m ade by a multiplying interferome- 
ter: DASI (Kovac et al]|2002 ). DASI had 13 single-mode anten- 
nas and performed pairwise correlation of signals across the Ka 
band, from 26 - 36 GHz. 

Several groups have studied the possibility of building a new 
generation of mm-wave interferometers spe cifically to sea rch for 
the small polarization signals in the CMB dTimbie et alJl2006t) . 
Compared to existing mm-wave interferometers, these new in- 
struments would have to do the following: 1) collect more modes 
of radiation from the sky by including more antennas (> 100 ); 
2) operate with broader bandwidth (~ 25%, corresponding to 
> 10 GHz), and 3) operate over a broader range of center fre- 
quencies, at least up to 90 GHz, to be able to detect and reject 



astrophysical foreground sources by their spectral signatures. 
Because pairwise correlation requires multiplying A^(A^ - l)/2 
signals together, conventional heterodyne interferometers face a 
significant challenge to increasing the number of antennas. 

On the other hand, adding interferometers have the ability 
to correlate large numbers of inputs over wide bands. Here we 
present an adding interferometer based on a Fizeau beam com- 
biner Co mbined with the optimal phase-shifting scheme de- 
scribed in lCharlassier et all (2009), this is a promising approach 
for measuring the B-mode polarization. The technique is com- 
patible with either coherent receivers (amplifiers) or incoherent 
detectors (bolometers). 

In a 2-element adding interferometer the electric field wave- 
fro nts from both antennas a re added and then squared in a detec- 
tor (i Rohlfs & Wilsonl2004l) . ( See Fig.[T]) The result is a constant 
term proportional to the intensity plus an interference term. The 
constant term is an offset that is removed by phase-modulating 
one of the signals. Phase-sensitive detection at the modulation 
frequency recovers both the in-phase and quadrature-phase inter- 
ference terms and reduces susceptibility to low-frequency drifts 
(1// noise) in the detector and readout electronics. The adding 
interferometer recovers the same visibilities as a multiplying in- 
terferometer: 



/(n, v)G(n)e'^'^-^Jiv)dvdn 



(1) 



where /(n, v) oc |£'o(n, v)p is the incident intensity and Eo(n, v) is 
the incident electric field as a function of position of the source 
on the sky, n, and frequency, v. G(n) is the primary beam (power) 
pattern of the antennas (assumed identical and for simplicity as- 
sumed independent of frequency), u is the vector between the 
centers of the antenna apertures, measured in wavelengths, and 
has components u and v. The frequency bandpass of the instru- 
ment is J(v). 

To combine signals from N > 2 antennas, we use a Fizea u 
beam combiner, a type of 'image plane' combiner dTraubl2000l) . 
This technique is analogous to the simplest interferometer in 1- 
d: the Young's double (or multiple)-slit interferometer (Fig.dHi. 
Schematics of Fizeau combination are shown in Fig. |5]). The 
beams from the apertures can directly illuminate an array of de- 
tectors, or, more typically, they pass through a lens or telescope 
first to reduce the size of the instrument. While Fizeau combin- 
ing is well-known, we stress here the fact that there are two 
types of path differences (and therefore phase differences) for 
rays traveling from a source to a detector: one path difference oc- 
curs outside the instrument, and the other inside the instrument. 
Compare this to a conventional interferometer (also shown in 1- 
d, though the extension to 2-d is straightforward) as shown in 
Fig.dUi, where rays only undergo a phase difference before they 
enter the antennas. (Note that long-baseline optical interferome- 
ters usually include a 'delay line' between the apertures and the 
beam combiner. This element introduces an equal path length for 
all rays entering an aperture. In contrast, the internal path differ- 
ences we are concerned with are different for each pixel in the 
focal plane.) 

Let us explore what this combination of path differences 
achieves. We start by noting that the 'external' phase differences, 
which are present in any interferometer, are the reason that the 
visibility function is a Fourier transform of the image on the sky. 
The visibility measured by a single baseline essentially selects 
one Fourier mode from the image. In the Fizeau system, we have 
an additional set of phase differences. Without loss of general- 
ity, we may assign a negative sign to the phases introduced inside 
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the instrument. Now, if we sum over both the phases, we get a 
Fourier transform followed by an inverse Fourier transform - but 
this is the image itself! Thus, Fizeau combination enables imag- 
ing in an interferometer. The image formed in the focal plane of 
the Fizeau combiner is equivalent to the 'dirty image' measured 
by conventional radio interferometers. This topic is discussed 
later in this paper in Appendix lAl 

In addition, we show that, given enough detectors on the fo- 
cal plane, we can extract some spectral information from the in- 
terference fringes and determine the visibilities in several sub- 
bands. The Fizeau system enables extraction of spectral infor- 
mation via geometry, without additional components like filters. 

Without spectral information an interferometer with a large 
bandwidth suffers from a large radial width of each pixel in the 
u-v plane. This is shown in Fig.©. Let v be the center frequency 
and Av the bandwidth. Then, a baseline of length B will measure 
the CMB power spectrum over a band centered on spherical har- 
monic number 



ttB nvB r-r 

{ - — = n Vu2 + v2 

A c 



where the width in ^-space is 



ttAvB 



(2) 



(3) 



As mentioned above, the additional spectral information that is 
available to us can be used to sub-divide the band in the u-v 
plane. We discuss this aspect in detail in ^ 

These advantages are studied in the context of two novel 
millimeter-wave interferom eters: MBI (the M illimeter-wave 
Bolometric Interferometer. (iTimbie et alJ l2006l) ^ and QUBIC 
dHamilton & Charlassieij|20 1oV In these svstems the apertures 
are replaced by an array of back-to-back single-mode horn 
antennas (see Fig.©). The inward-facing horns illuminate a 
Cassegrain telescope that combines the beams at the focal (im- 
age) plane, which is tiled by an aiTay of bolometers whose di- 
mensions are small compared to the overall dimensions of the fo- 
cal plane. Between the outward and inward facing horns are elec- 
tronic phase modulators which can be operated independently. 
They are used to modulate the interference fringes that appear 
on the focal plane in such a way that fringes ca used by different 
basel i nes can b e distinguished from each other (Charlassi er et alJ 
l2009tlHvland et al. 2009) . In principle this technique could work 
with multimode horns or antennas as well. 



3. The Fizeau combiner output and its relation to 
visibilities 

In this section, we study the output of the adding interferometer 
and its dependence on instrument parameters: number of detec- 
tors on the focal plane, number of antennas, etc. and relate this 
output to visibility from an interferometer, in order to describe 
bandwidth splitting in the following section. In what follows, we 
denote the output at the detectors as O. A simple adding interfer- 
ometer with two antennas/apertures (single baseline) is shown in 
Fig O- A generalized Fizeau adding interferometer is shown in 
Fig.®. 



3.1. Simple Interferometer 

In an adding interferometer, electric fields are added and then 
squared at the detector. If Ei(n, v) and E2(n, v) are electric fields 



Primary Beam 
G(n) 




Square-Law 
Detector 



G(n)[Eo^+Eo^cos((p+A(p(t))] 



Fig. 1. Schematic of an adding interferometer with N = 2 antennas. 
At antenna A2 the electric field is Eg, and at A[ it is iioexpi^, where 
if = 2;ru-n and |u| = \B\/A. |B| is the length of the baseline, and (B-n)/|B| 
is the angle of the source with respect to the symmetry axis of the base- 
line. (For simplicity consider only one wavelength. A, and ignore time 
dependent factors.) In a multiplying interferometer the in-phase output 
of the correlator is proportional to cos tp. For the adding interferom- 
eter, the output is proportional to E^ + E^ cos((^ + A(p{t)). Modulation 
of A(f(t) allows the recovery of the interference term, Eg cos (/>, which is 
proportional to the visibility of the baseline. 



Horn antennas 



phase 
modulators 

inward- 
facing antenna: 



, skipping (N-2) antennas . . , 
Secondary mirror 



Primary Beam 
G{n) 



G{xl 




Primary mirror 



Focal plane detector array 



Fig. 2. Block diagram of an adding interferometer with N > 2. Each 
phase shifter is modulated in a sequence that allows recovery of the 
interference terms (visibilities) by phase-sensitive detection at the de- 
tectors. The signals are mixed in the beam combiner and detected. The 
beam combiner can be implemented either using guided waves (e.g. in 
a Butler combiner) or quasioptically (Fizeau combiner), as above. The 
top triangles represent corrugated conical horn antennas. For the case 
of an interferometer using coherent receivers, amplifiers and/or mixers 
could be placed before the beam combiner. 



incident at antennas 1 and 2 respectively, then, denoting the pri- 
mary beam of the outward-facing antennas by G(n), an adding 
interferometer will detect 



G(n)|Ei(n,v) + E2(n,v)|2 



(4) 



from a certain direction n. Now Ei(n, v) and E2(n, v) differ only 
by a phase factor tp = 27rB ■ fi/A, so that we can write the electric 
fields as E()(n, v) and Eo(n, v) exp(/^), and the output at a single 
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frequency for the simple arrangement shown in Fig.© is 



^4l] andl4!2l For now, Eq.® reads 



G(n)|Eo(n, v) + Eo(n, v) exp(/^)|2 
G(n)|Eo(n,v)|2(l +cos^) 
G(n)|Eo(n),y|2(l+!R[expOV)]) 



(5) 
(6) 



ZSqCA, v) oc /(n, v) where /(n, v), the intensity, is a linear combi- 
nation of Stokes parameters (see Lin & Wandelt 2006, e.g.)- To 
isolate the interference term, a time-varying phase modulation 
is applied, as shown in Fig.©. Then, the demodulated output is 
given by 



G (n) /(n, v) cos 1^ = G (n) /(n, v)^[exp(!V)] (7) 



where - IttB ■ h/A - Inn ■ n. 



3.2. Fizeau Interferometer 

The general arrangement of a Fizeau interferometer is shown in 
Fig-©- Let Ba be the baseline formed by antennas at and 
A,, such that = A^ - A; (k e [I . ..N(N - l)/2] and p, I e 
[1 . . .A^]), and n is a direction in the sky, as shown in Fig.©. 
Then, Uk = B^/zl, so that the external path difference is (Ap - 
A/) ■ n = Bjt ■ n, as shown in Fig. Q. 

The difference between this arrangement and that of the sim- 
ple adding interferometer in Fig.Q is the set of internal phase 
differences introduced inside the instrument due to the fact that 
there is more than one detector, and the geometry of the arrange- 
ment. Let Xj denote the position of the jth detector on the focal 
plane. Then, the internal path difference between the rays from 
antennas p and I, (which form baseline k), is ^^^(x) = ■ B^ as 
shown in Fig.®. Thus, the intensity at the detector at x^ con- 
tributed by baseline k is 



2;r 

G (n) /(n, v)%[exp i{ipk + —Xjt)]. 

A 



where (fk - 2jtBi! ■ h/A - 2n\ik ■ n. 

If we denote the output at the jth detector from the A:th base- 
line as Ojk, and integrate over all directions in the sky and the 
bandwidth Av, we get 



Ojk^ j j G(n)/(n,v)!R[exp!(^, 



27r 

'k + —Xjk)]J(v)dndv, 



exp(i—xjk)x 
Aq 



J Jg(^) I(n,v)exp[iiipk)] J{v)dvdh 



(10) 



[ 2;r 

%\e.x-p{i—Xjk)x 
Aq 



J J G(n)/(n, 



v) exp [/27rUi: ■ n] J(v)dvdn 



(11) 



The quantity indicated by the underbrace in eq.dTTTi is the 
Visibility (defined in eq.(l)) from the A:th baseline, V(Uk) = Vk. 
We can now denote the 'internal phase differences' (2n/A)xjk as 
(pjk, so that 

Ojk^n[sxp(i<Pjk)Vk]. (12) 

In AppendixlAl we show that by integrating over the entire focal 
plane, we can recover the image convolved with a "dirty beam", 
as in radio interferometry. 

In the next section, we consider the net signal from a single 
baseline and describe how it can be split into 'sub-bands.' For 
the sake of simplicity, the index k corresponding to baseline k is 
dropped, and all equations in ^hold for every baseline. 

4. Spectral information from an interferometer 
using a Fizeau approach 

4.1. Preliminaries 

The output measured at the detectors in a Fizeau interferome- 
ter contains the following phase information integrated over the 
entire bandwidth: 



1 . phase introduced because of the path difference between any 
two rays that arrive from the same part of the sky on the two 
outward-facing antennas that make up a baseline; and 

2. phase introduced because of the path difference between any 

(8) two rays that arrive from two different antennas on to the 
same point in the focal plane. 

The phase in point[T]is due to the fact that we are considering a 
radio interferometer, and so the visibility that we measure must, 
by definition, include this phase. However, the phase in |2] above 
introduced by the beam combiner needs to be factored out to 
recover visibility from each bolometer. If there were a way to 
calculate the 'internal phase' introduced by the beam combiner 
over the whole bandwidth, then all we would need to do is to 

(9) divide the output at each point in the focal plane by this 'internal 
phase', and we would get visibility directljQ. 



where J(v) is a bandpass weighting. There is an integration over 
detector area as well in eq.®, which will be implicit until 34.31 
where we explore the effect of the detector area and make a 
suitable and practical approximation . Notice that xjk does not 
depend on the direction n. For now, we make the (crude) ap- 
proximation that A = Aq (the central wavelength), so that this 
'internal phase factor' may be taken out of the integral over v 
as well. As discussed in 34. II and 14.21 this crude approximation 
need not be made. We will instead choose to view the visibilities 
as 'averages' over a certain 'sub-bandwidth', in which case this 
'internal phase' factor is an average over this 'sub-bandwidth'. 
We will discuss a more efficient way to deal with this issue in 



4.2. 'Sub-band Splitting' 

We can think of the effect of the instrument on the visibilities 
in the following way. Let us divide the entire bandwidth of the 

' A range of values of v will produce a range of f's, or a band in 
^-space (as discussed in ^ see eqs.l i2l3t ). A finite-bandwidth interfer- 
ometer thus measures what is called a 'bandpower' instead of a single 
value of the power spectrum at one value of £. But the power spec- 
trum is jus t the variance of the visibilities for a circle (ring) in the u-v 
plane (see iMal j[2007l : ICharlassier et alj|201(]|. Fig. 5.8, Figs. 1&2 re- 
spectively). And so we get different bandpowers for the same baseline 
and orientation but for different frequencies. 
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instrument into m sub-bands and let vi, V2...v,n be the centre- 
frequencies of each one. Then for one baseline, one orientation, 
and one detector position, these will correspond to visibilities 
Vi, V2---Vm and to phase differences (pji, (f)j2...4>j,„ (where j rep- 
resents the detector). If we represent the output at the jth detector 
as Oj then we get (as in the previous section, but dropping the 
index k for baseline, since we're considering just one baseline): 

m 

O, = 2!R[y„exp/0,„]. (13) 

a=l 

Given just one detector, it is impossible to extract every V„ for 
every sub-band, even though we know precisely what the (pja's 
are. However, if we have m detectors, then we can easily write 
the following system of equations for each baseline: 

01 = %[Vi exp i^u + V2 exp /0i2 + ■■■ + ¥,„ exp /(^i„,] 

02 = Vi exp i^2i + V2 exp /022 + ■■■ + V,„ exp i^2m] 

Om = tR[yi exp /0„,i + V2 exp /0„,2 + ■ ■ ■ + Vm exp /0„„„]. (14) 

Now, as discussed in the preceding section and in Fig.(l), we 
can apply a unique phase shift to each baseline. If we denote 
this phase shift by Aip (it is understood that = A^(f)) and the 
output after applying the phase shift as O'y . . .0',„, then we get 
another set of m equations: 

0\ = !R[yi exp i{(pi 1 + A(p) + ■■■ + ¥,„ exp i{(pi„ + Aip)] 
O2 = 'KlVi exp /(021 +A(p) + --- + V„, exp /(02m + A^)] 

0'„, = '^[Vi exp/(0,„i +A<p) + --- + V,„ exp/(0„„„ + A^)]. (15) 

This is a system of 2m equations with 2m unknowns - 1R[Vi], 
%[V2] . . . %[V,„] ...3[Vi], 3[V2] . . . 3[y,„], and so we can solve 
for the values for each one of these 'sub-band visibilities'. The 
beam combiner thus achieves far more than just separating the 
real and imaginary parts of visibilities. 



4.3. Effect of finite detector size 



We mentioned in ^4.1l that if the 'internal phase' is known, sub- 
band visibilities may be computed/estimated. However, calcu- 
lating this internal phase is not easy, since integration over the 
bandwidth complicates the calculation, as seen in ^3.21 eq. (fTTl i. 
where the 'internal phase' factor was taken out of the inte- 
gral to define a visib ility using a crude approximation (see 
ICharlassier et al. 2010. §2.2). The internal phase differences in 
the Fizeau interferometer create fringe patterns on the focal 
plane. This is exactly the same as saying that the visibilities 
in each sub-band are modulated by a fringe which depends on 
baseline length. To extract these visibilities, we need to separate 
the fringes. In order to do so, we need to realize that what we ob- 
serve at every detector is the visibility on the sky times the fringe 
summed over the area of the detector as well as bandwidth. 

The fringe pattern is different for every frequency in the 
bandwidth. Visibility is also different for different frequencies. 
Since we can compute the fringe corresponding to each fre- 
quency, it is also straightforward to sum up these fringes over 
a small 'sub-band' over the area of a single detector The output 
at each detector from a baseline is known, and this can be written 
as a sum over the product of visibility for a 'sub-band' and the 
fringe for the respective 'sub-band'. This system of equations 
can be solved for each baseline to yield the sub-band visibilities. 



We have demonstrated this in ^4.21 Effectively, this amounts to 
not using the crude approximation in ^3.2l eq. (fTTl l and instead 
writing the output as a product of a sub-band visibility and a 
fringe. 

In the discussion so far, we have assumed that the collecting 
area of each detector is negligible, and we completely ignored 
the effect of the fringe pattern. Let us account for these effects in 
the following way. Let A be the effective collecting area of each 
detector Let / (x, Vq.) be the value of the fringe pattern (see dis- 
cussion at the beginning of this section) at a point on the focal 
plane x and in a frequency sub-band marked by a. Then, equa- 
tions (114115b become 



III n 

Oj I y„expi0;a(x)/(x, v„)<i^x 

"' r r 

O'j = J]^ I Va exp i(4>j„ (x) + Aip)f (x, V J d^x 



(16) 



where it is understood that integration is done over the area of 
the detector 

This leaves us with an issue - that of deconvolving the Vs 
from the integrals. However, if the area of the detector is small 
compared to the width of fringes, then we can assume that the 
phase differences remain roughly constant over the collecting 
area of one detector, so that we may write 



(x)F(x, 

a=l 
m 

Oj = A ^ 51 [y, exp /(0,„ (x) + Aip)F (x, y,)] (17) 

where F (x, v^) represents an "average" value of the fringe pat- 
tern. 

Equations ( fTTl i again have 2m variables and can be solved to 
get 2m quantities: the real and imaginary parts of m visibilities 
over the bandwidth. 

Application to CMB cosmology: The QUBIC collabora- 
tion is im plemen ting the technique described in this paper 



tion is im plemen tmg ti 
(ICharlassier et al.,i2010h . 



5. Conclusions 

1 . The Fizeau system makes it possible to recover spectral in- 
formation without the need for filters. 

2. The Fizeau system acts naturally as an imager 

In addition, by introducing phase modul ators discussed in 
(iHvland et al.l l2009t Iciiarlassier et al.l l2009h . we can measure 
visibilities for all baselines in a Fizeau system. 

While it is possible to divide the bandwidth into many differ- 
ent sub-bandwidths, it isn't possible to do this indefinitely. The 
beam for a single antenna determines the FOV of the instrument 
and limits the resolution in the u-v plane, as shown in Fig.(l3]l. 

It is also possible to operate the interferometer simultane- 
ously as an imager The additional modulation mentioned above 
opens up a range of possibilities, including the simultaneous 
measurement of visibilities and images. This is described in 
Appendix lAl 

In conclusion, the Fizeau system introduced here is poten- 
tially powerful tool for astrophysics: it could allow the recov- 
ery of more information than is possible with traditional inter- 
ferometers or imagers and does not need significantly more re- 
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sources to build. Its application in CMB cosmology is straight- 
forward and can be demonstrated in future version of QUBIC 
(Hamilton & Charlassier 2010). 

Acknowledgements. We thank the membens of the MBI and QUBIC collabora- 
tions for many fruitful discussions on bolometric interferometry. 
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Sub-band division of pixel widtli=FOV size 



V 




U 



Fig. 3. The u-v plane coverage of one baseline of an interferom- 
eter for a single pointing in a single baseline orientation angle. 
Radial spread in a single pixel in the u-v plane due to bandwidth 
is shown. Resolution in the u-v plane is determined by the pri- 
mary beam or size of field-of-view. The minimum size of each 
sub-band is also determined by this resolution, as shown. 
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Fig. 4. A Fiz eau adding interferometer used as a beam combiner for long-baseline optical interferometry (figure reproduced from 
iTrau b (2000)). The baseline B12 is the separation of the apertures, which can be widely spaced, and determines the angular resolution 
of the instrument. The Fizeau beam combiner has a different and much smaller baseline length (not labeled). The 'external' phase 
differences are marked B12 ■ n. 




Fig. 5. A simple 1-d Fizeau system, similar to a figure in lTraubl(l2000l) . O denotes the origin of the co-ordinate system, and is the 
center of the observation plane. x,'s denote positions of detectors on the focal plane, and A/'s denote positions of apertures, where 
/ G [ 1 . . . A^] . Notice that there are two sets of phase differences, marked (Ap - A/) ■ n = ■ n and Xjk{x) - Xj ■ (Ap - A;) for external 
and internal phase differences respectively. The geometrical variation of Xjt(x) = ■ (Ap - A/) allows "sub-band splitting". 
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Appendix A: Imaging in a Fizeau system 

In eq.(|9]l, if we now integrate over the focal-plane area, omitting 
the integral over bandwidth, which is implicit (and remembering 
that if = 27rB ■//!): 



O 



G(n)/(n,v)exp 



i2;7r- 



B ■ n - jc, 



dnd'-iL, (A.l) 



where we have changed the sign on the 'internal' phase differ- 
ences. This can be done without loss of generality, since the 
internal and external phase differences are independent of each 
other 

Let us consider just one term in the expression (...): 



O 



v) exp 



i2;r- 



B n 



c/nexp 



-Yin— 



d^x (A.2) 



If we include the effect of the primary beams of the inward- 
facing antennas (G(x)) and adopt / = /(n, v). 



O 



j G(x) J G(n 



)/exp 



i27r- 



B n 



dh exp 



-i2;r^i 



d^x 



(A.3) 



The quantity in underbrace is clearly a fourier transform, and the 
expression can be written as 



O 



I 



G(x)g(G/)exp 



-i27r- 



d^x 



(A.4) 



If the distance from the inward-facing antennas to the focal plane 
» the collecting area for each bolometer. 



(9 = r'(G5(G/)) 



(A.5) 



The beam needs to be deconvolved from the above expression in 
order to obtain an image from the instrument. 

Now, eq (IA.4l i can be split up over the focal plane: 



O 



N „ 



Gg (G/) exp 



-i2;r^ 



d^x 



(A.6) 



thus tells us that the dirty beam for the image produced by the 
Fizeau combiner is more involved than the traditional interfer- 
ometer dirty beam, but remains conceptually equivalent. 

There are two assumptions inherent in the foregoing discus- 
sion: 

1 . The focal plane is large enough to receive most of the power 
from the inward-facing antennas 

2. There are no "blank" areas on the focal plane for which the 
incident power is not absorbed by a bolometer 

It is possible to detect the correlated signal from a pair of an- 
tennas as well. In order to separate every unique baseline, time- 
varying phase-shifts can be applied to the signal at the base of 
the skyward-facing antennas. The "coiTelated signal" from each 
pair of antennas is simply the visibility from a baseline, with one 
crucial difference: there is an "internal phase" added to every vis- 
ibility due to the relative positions of antennas on the observation 
plane and detectors on the focal plane. These phase differences 
are geometrical. Since the output from all antennas is inci- 
dent on every detector, we can say that the output from a single 
detector contains information about A^(A^ - l)/2 visibilities. 



where 1 . . . are labels for bolometers on the focal plane. 

Each of the bolometer outputs then represents a pixel in im- 
age space. The total number of pixels depends on the resolution 
of the instrument, and not the number of bolometers on the focal 
plane. Therefore, if the number of bolometers on the focal plane 
are greater than the number of pixels in the image, we need to 
"repixelize" the image obtained, so that all pixels are indepen- 
dent of each other 

In general, this is how the beam is convolved with the image 
on the sky for the Fizeau beam combiner: 



O = r'(G5(G/)) 

= [(r'G)*(G/)] 

In traditional interferometry, eq. iA.li would read 
O = r'(5(G/)) 

= (gG) * m 

= GI 



(A.7) 
(A.8) 



(A.9) 
(A. 10) 
(A. 11) 



In eq. dA.lOb . the u-v space beam gG needs to be multi- 
plied by u-v coverage, which is a "mask", say W(m, v). Then, 
5"' (gG X 9J!(m, v)) is called the "dirty beam" in traditional in- 
terferometry. In eq. (IA.7l i. the factor 9[R(m, v) is included. Eq. dA.SI l 



